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APPLICATION  OF  MIXED  STRATEGIES  FOR 
IMPROVED  MISSILE  GUIDANCE 

Grant  No.  AFOSR  89-0528 

SUMMARY 

The  main  task  in  the  reported  period  nas  been  the  preparation  of  the 
for  a  non-dimensional  sensitivity  analysis.  This  sensitivity  analysis  will 
to  generalize  the  results  obtained  so  far  in  previous  phases  of  our  resear: 
to  perform  a  parametric  investigation. 
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INTRODUCTION 
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The  work  on  the  task  of  non-dimensional  sensitivity  analysis  started  in 
November  1990  in  parallel  with  the  preparation  of  the  report  on  the  last  year’s 
activities.  The  progress  of  the  work  has  been  slower  than  expected  for  two 
separate  reasons:  the  need  to  introduce  and  train  new  research  personnel  and  the 
perturbations  caused  by  the  Gulf  War. 

At  present,  we  are  in  the  last  phases  of  preparing  and  validating  the 
computational  tools  for  the  parametric  investigation. 

This  Interim  Report  describes  two  topics: 

a.  The  principles  of  the  planned  non-dimensional  sensitivity  analysis. 

b.  The  practical  steps  taken  in  preparing  and  validating  the  computational 
tools  for  the  practical  investigation. 

NON-DIMENSIONAL  ANALYSIS 

In  order  to  generalize  the  numerical  results  of  any  study,  it  is  a  common 
practice  to  express  them  in  a  proper  non-dimensional  form.  In  order  to  do  so  the 
following  steps  have  to  be  taken:  Scaling  of  all  variables  of  the  problem,  based 
on  a  judicious  selection  of  scaling  parameters  of  some  physical  significance. 
Definition  of  the  "similarity  parameters"  of  the  problem. 

The  present  investigation  effort  has  been  guided  by  some  well-known  results 
of  "linearized  missile  guidance"  analysis  and  their  extension  as  reported  in 
Refs.  1  and  2. 

Since  the  study  involves  point-mass  models  only  the  last  two  of  the  three 

basic  dimensions  (mass,  distance  and  time)  have  to  be  scaled.  In  our  dynamic 

model  the  missile  is  assumed  to  have  a  first-order  transfer  function  with  the 

time-constant  x  .  The  value  of  this  time-constant  is  selected  therefore  as  the 
P 

unit  of  "normalized"  time. 


2 


t  =  t/(t) 


ref 


(1) 


(t)  _  =  T 

i  et  p 


(2) 


The  unit  of  "normalized"  distance  was  selected,  based  on  the  fact  that  the 

study  is  concerned  with  maneuvering  targets,  as  a  product  for  maximum  lateral 

2 

acceleration  of  the  evader  (a_)  and  t  .  (In  the  "classical"  linear  guidance 

£  max  p 

analysis,  miss  distances  due  to  target  maneuver  are  proportional  to  this 
factor.  ) 


x  =  x/x 


(ref ) 


,3) 


(x)  =  t  (a  )  (4 ; 

ref  p  E  max 

Using  the  above  defined  scaling  transformation  keeps  the  following  non- 
dimensional  variables  unchanged: 
a.  The  pursuer-evader  speed  ratio  {ft) 


V  V 
3  A  2  =  _P 

VE  V 

E 


(5) 


b.  The  pursuer-evader  maximum  maneuver  ratio  (/a) 


a  (aD)  ( a  _, ) 

A  P  max  _  P  max 

M  "  TIT)  Tall 

E  max  E  max 


(6) 


c. 


The  normalized  time  of  flight  (9^.) 


0.  £ 


o  1 


f  (  t  )  v  T  - 

ref  c  p  V 


(7) 


d.  The  "dynamic  similarity  parameter"  (Refs.  1,2)  defined  by 


.  (VI  (a  ) 

.A  E  max  E  max 

y  =  — v -  TP  =  — - 

c  r  \r 


(8: 


Using  the  time  scale  change  given  by  (1)  all  frequencies  are  normalized 
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(9) 


0)  =  t  ~  )  =  WT 

ref  p 

This  scaling  equally  refers  to  the  actually  assumed  periodical  random  phase 
maneuvers  and  electronic  jinking  as  well  as  to  the  shaping  filters  modelling 
these  stochastic  processes.  The  spectral  density  of  the  white  noise,  used  as  the 
input  of  the  shaping  filter  has  to  be  also  rescaled  according  to  the  units  of 
the  considered  disturbance. 

In  our  previous  investigation  three  types  of  disturbances  were  considered: 

1.  The  "glint"  noise  of  the  radar  receiver,  modelled  as  measurement  noise. 

2.  The  periodical  random  phase  target  maneuvers. 

3.  The  "electronic  jinking". 

The  last  two  were  considered  as  "process  noises".  As  a  consequence,  another 
important  similarity  parameter  has  to  be  identified:  the  ratio  of  the  two 
process  noises,  (in  other  words  the  "normalized  jinking"  parameter) 


(W)  =  (U)  /x2.  . 

max  max  p(a„J 

E  max 


(10) 


Moreover,  the  "glint",  which  is  proportional  to  the  physical  size  of  the  target 
aircraft,  must  be  scaled  as  all  the  other  variables.  Since  the  dimension  of  the 


glint  spectral  density  $  ^  is  [m  sec],  the  normalized  spectral  density  s 


houid  be 


=  $  /x  (a  ) 

gl  p  E  max 


(11) 


Keeping  all  the  above  defined  similarity  parameters  constant  should  lead, 
according  to  the  theory,  to  the  same  normalized  miss  distance 


M  =  R( t  )/x2(a  ) 

f  p  E  max 


(12) 


against  the  same  evader  strategy. 

In  order  to  guarantee  the  same  operational  efficiency,  i.e.  the  same  single 
shot  kill  probability,  the  lethality  range  of  the  missile’s  warhead  has  to  be 
kept  proportional  to  the  expected  miss  distance.  This  can  be  done  if  the  non- 
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dimensional  warhead  lethality  parameter,  defined  as 


*1 


R?  =  R/r^a  ) 

i  t  p  L  max 


(13) 


is  kept  constant. 

In  summary,  we  have  a  set  of  seven  non-dimensional  similarity  parameters 

13,  p,  y* ,  0.,  (W)  ,  d>  ,  ,  Rt,  as  a  basis  of  a  parametric  sensitivitv  analysis. 

f  max  gl  l  J 

In  theory,  for  the  same  set  of  these  parameters  similar  non-dimensional  results, 
expressing  a  large  family  of  dimensional  results,  arc  to  be  expected. 

The  number  of  parameters  being  quite  large  it  seems  to  be  useful  to  keep 
some  of  them  constant  and  to  carry  out  the  parametric  sensitivity  analysis  by 
varying  only  a  few. 

For  the  planned  sensitivity  analysis  the  following  parameters  are  going  to 

be  kept  constant: 

the  speed  ratio  fi  =  V_/V 

P  e 


the  dynamic  similarity  parameter  y* 


x  (a  ) 
P  E 


max 


Ro  1 

the  normalized  flight  time  — 

c  p 


The  remaining  parameters  p,  (W)  ,  d>  ,  ,  R*  belong  to  two  different  groups.  For 

max  gl  i  or 

a  given  scaling  transformation,  based  on  x  and  (a_)  ,  p  and  R*  can  be 

p  E  max  i 

considered  as  design  parameters  of  the  missile,  while  (W)  and  d>  ,  are  related 

max  gl 

to  the  evading  aircraft.  The  last  two  both  depend  on  the  physical  dimension 

of  the  aircraft  and  have  to  be  varied  together  accordingly.  They  signify  the 

relationship  between  disturbance  created  by  the  evader’s  maneuverability  and  the 
disturbance  due  to  ECM  or  glint  respectively. 

Another  important  distinction  between  the  two  pairs  of  these  parameters  is 
that  the  first  two  (p  and  Rj| )  do  not  interact  with  the  estimation  process  and 

therefore  no  change  in  the  "optimal  estimation"  is  expected  with  their 

variations.  However,  if  the  ratio  of  the  different  disturbances  is  changed  it 
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most:  certainly  will  require  a  modification  of  the  estimator. 

In  this  way  the  sensitivity  analysis  will  be  separated  in  tow  parts: 

a.  sensitivity  of  the  guidance  performance, 

b.  sensitivity  of  the  estimation  quality. 

This  separation  will  also  allow  to  obtain  the  correlation  between  the 
estimation  error  and  the  miss  distance. 

IMPLEMENTATION  ISSUES 

There  are  two  options  to  carry  out  the  above  outlined  nr.r.-d  i  u-r.'i  1 

parametric  investigation  The  first  one  is  to  use  the  existing  ccT.put'-r  cod--* 
which  is  written  in  dimensional  form,  by  carefully  adapting  the  input  in  a  v;. 
that  the  required  similarity  is  preserved.  The  second  option  is  to  rewrite  the 
computer  code  in  non-dimensional  variables. 

In  both  cases  careful  validation  tests  has  to  be  carried  cut  in  order  t 

gain  confidence  in  the  method  of  the  analysis  and  the  results  Validation  for 

each  case  is  via  comparison  to  some  existing  dimensional  data. 

Presently  we  are  engaged  in  work  using  the  first  option  and  have  already 
discovered  already  several  problems  to  be  addressed.  In  parallel  with  tne  above- 
effort  a  new  non-dimensional  computer  coae  Is  being  generated 

CONCLUDING  REMARKS 

In  this  interim  report  the  analytical  basis  of  a  non-dimensional  parametric 
sensitivity  analysis  is  outlined  and  the  two  options  to  carrv  out  the  task  are 
indicated.  Presently,  work  is  being  done  parallel  for  both  options. 
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